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INTRODUCTION 


The purpose of this research was to investigate some 
of the effects of boundary wall perturbations on finite 
amplitude standing waves іп а rigld-walled rectangular 
cavity.The investigation was prompted by an examination 
of the experimental results of Coppens and Sanders[3], 
the research of DeVall[5] and of Kilmer[4] ,which suggested 
the existence of the excitation of modes other than 


those belonging to the family of the driven mode. 





2.BACKGROUND 


A plane elastic wave travelling in a non-dissipative 
fluid will change waveform as predicted by the relevant 
hydrodynamic equations [8] ,[7] ..ır the problem is extended 
to absorptive media,only waves of relatively high amplitude 
will change waveforms appreciably. 

TE ала Россосаацаве sehool,Coppens andySanders 
(3] ,Kiimer [ 4] , and DeVall [5] have dealt with the study 
of finite amplitude waves in rigid-walled rectangular cavities. 

One interesting result of these cavity experiments was 
the appearence of excitations of modes which were not 
family members of the driven mode.For example,assume a rigid 
cavity of dimensions ILE is driven acoustically at 
frequency w,the resultant pressure standing wave is of the 
тот 

cosk x cosk y cosk,z cos(wt+8) oca) 
where k zNT/L. s k AN а LA 2) 
and N,M,L are integers.Eq.(2.1) can be represented by the 
notational shorthand 

ШІЛ Б/у,8) 
If the cavity is driven to excite the (0,M,0) mode,then the 
family of standing waves consist of all of those of the form 
(0,nM,0/nw,8. ) when DRM NE 

As it is stated in [3),"The standing waves which can be 
excited in any real cavity deviate from the predictions of 


the linear wave equation with ideal boundary conditions 


1:0 





"for the following reasons: 
(a).The presence of boundary-layer losses at the cavity 
surfaces yields a dispersive contribution to the wave equation. 
(b).Geometrical irregularities alter the effective 
dimensions of the cavity. 
Both of these mechanism can be treated as equivalent as 
long as the shift in frequency are so small that the actual 
resonances are close to the theoretical values resulting 
from the classical model." These are treated by assuming 
the dimensions are exact,and the apparent phase speed is 
determined on that basis. 

The resonance frequency for each standing wave is 
defined as [ 3] 


25612 
kd (2.3) 


where k's are given by Eq.(2.2) and nis a shorthand for 


l 2 2 
DEC > (nj k^ * (n, 


the set (By ony nz) where n 


the apparent phase speed appropriate for that frequency. 


non 


go My n; are integers ‚and en 1s 


The non-linear wave equation applicable to this 
problem can be obtained as follows: 
The continuity equation for wave propagation in Eulerian 


coondinates is 


E 
mer) + e (2.4) 


БЕШ ШЕ ОО sincan be written im terms of the condensation, 
s=(/-A)/A as 
> 95 
V.E(1+5)4] +35 =0 (2.5) 


The equation of motion in Eulerian coordinates for a 


contar ned* iscous fluid is 


ДЕСІ 





Pee a v(v.ü)-7vxvxü TODAT 


Qu OD 
where 
2. 
Pa 4i (L) = 46 [1099 UE s] (2.7) 


ODAT=Other dispersive and absorptive terms arising from 
boundary effects. 
Eq.€2.6) can be rearrange in the form of 
A у) ++ vp=+ fu De® 

E. the operator A describes those physical processes 
leading to absorption and dispersion. 

One can write Vxl-o and therefore U-vö „where Û is the 
Cl ОООО СП лд based оп Еле arrotational velocity 
assumption.Hence, fus £L VQ „Кер eng u=VO and using the 


Gemdemsation,b£q.€2.8) cam be written as 


Уф ++ Tv) + fe La 9 (145) E = 720 eo 


I sEcheschelp of Eq. 0225)9and(2.9), and after a good 
deal of manipulation the non-linear wave equation may be 


approximated in terms of EA. potential 


ep a, $ = [(và)* Ez Ed. 2] 





on.) 
where 
шй. = СП + 2-5 
It should be noted that if the fluid is lossless and 
-Cg then (2.10) reduces to a previously known non-linear 


wave equation [ 9] 


12 





|- 
exe 
NS 
et o: 
bs 
و ا‎ 


2 ° = 
c 0^6 s 4 [(70ф) + = Cz 
(2.104) 

In order to express the approximate non-linear wave 
equation in terms of acoustic pressure and particle velocity, 
one can rearrange the Eq.(€2.9) in terms of p and Ù and 
mS Chat eJUaltlolh ШЕШ (2.10).Тһпе result 15 a quadratically 


non-linear wave equation [2] 
| ú x 
cO GE.) 9 - T gu [* (25 2а) +(e) J 


272 P \2 ú y? 
те” Ї (Z) - (ë) ] 


£ 11) 


2 a 
If it chances to be that ( aie | and (=) nearly satisfy 
fo Ce Co 
the wave equation, СУП Е | then on the RHS of Eq. 
2 
. 11) ES E = апа (2.11) becomes 
2 y P En и 
le HF (ZY+ (=) 1 om) 
Pines ITE I happens that Ko : 
PP 2t (Z) = = iR E 


as is true for solutions to the wave equation separated in 


ЕЕ Ешан соозапатеє "then } 


cp oL £c] ie [ = D (Z= C y] 2339 
EE ts. Ts crus om hor ear tesian coordinates.) 
Ip f ic i bh Joi Or Eq.(2.13) is the 
classical,linear wave equatton pertinent to the system 
ене Rio Can ве interpreted as a forcing function 
consisting of a three-dimensional spatial distribution of 


phase-cohorent sources.In a second-order perturbation theory, 


> 





МЕЕ f cano аса 1S Obtained from the classical (first- 
order) solution of the acoustic problem.The second-order 
perturbation solution describes the nonlinearities resulting 
from the self interaction of the classical solution.Higher- 
order perturbation solutions consider the interaction of the 
fom—“inear solution with itself,and the forcing function is 
КО осе Of products of Dorn classical and nonlinearly 
весе terms. 

Thus ,if a system is driven at frequency w,the nonlinear 
Pn "eguation 029079 9 will force the existence of 
all integer multiples nw of the driving frequency and the 
Eee on Ton must contra Ll harmonics of the input 
frequency.In a closed cavity,each of those nonlinearly 
generated waves whose frequency lies close to the resonance 
frequency of a standing wave of the cavity and whose 
associated spatial function matches that of the standing 
Wave can be strongly excited.Just how strongly will depend 
on the quality factor Q for the particular resonance and 
the difference between the resonance frequency of the 
standing wave and the harmonic nw..... 

"Consider two limiting cases. 

ШІРЕНЕ ЕЕ Foreinez2-merion does not have its frequency 
nw close to w ‚this standing wave is being forced at a 
frequency far removed from its resonance.Since this yields 
the inequality > 

D 
| Co OP | » It BP | (2.14) 


Eu bc ismored in....Eq.(2.13). 
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Сс nwe,w ,then the standing wave is being forced 
near resonance,and losses must be retained in..."Eq.(2.13). 
"The value of С. сап be determined from the apparent dimensions 
of the cavity and the measured resonance frequency Wo 

"The losses are described by the measured 2 of the 
resonance.This means that the linear-wave equation operator 
for the P. can Бей ығып ^. г 

202 nw 

CoD * Ж $5 a с, Ее Qa ?t 16) 

"Comparison. of cases (1) and (2) reveals that the 
responce of the cavity when nw w is order of 1⁄ e area 
to that when nwa w -Thus for the high-Q resonances usually 
СО ОШ ered IM Gavities Еп гісіс walls,the components in 
the forcing function which excite standing waves far from 
resonance can be ignored compared to those components 
exciting standing waves near resonance.... 

"The non-linear,coupled,transcendental equation appliciable 


to this problem can be expressed as 


COS 
x [ois] (0,9, 0 =N n Bo, cose [1/2 > Pia 2. +Z _ 


n 
cos 
dst Т. ЖҰЗ ТТЫ?) 


for all n51l.The values of ШЕ апа Wo must be determined from 
the infinitesimal-amplitude behavior of the cavity.The Mach 
number Mo and driving frequency w are known and Ng has the value 


N.=1/2 for a one-dimensional standing wave 
1/4 for a two-dimensional standing wave 


1/8 for a three-dimensional standing wave." 


25 





Ra И ШЕ ШЕ сел СО Cent Of nth harmonic component, 
normalized such that R z1. P. is the phase angle of the 
nth harmonic component ‚where f, 20. and Fe se ele ibe or 


is given by[3] 





canoe >t om 
where A EU ON neun) for Xn ¿1 Сао 
А м 
n 


Equation (2.16) can be solved by a method of successive 
approximations on a digital  computer.This has been done by 
[3] and[5] and both decided that the theoretical model can 
be used to identify the modes of a non-ideal,rigid-walled 
cavity provided quantities e, to be defined later are 
sufficiently small.The theoretical model in its present 
Bormardılesceor account for спе excitation of modes other 
than those belonging to the family of the driven mode. 

This excitation was observed to occur only in the case of 
nearly degenerate modes.It is believed to be caused by 
some linear coupling mechanism within the cavity. 

The purpose of this research is to see if the presence 
of wall irregularities can explain how non-family members 
may be strongly excited,and to present an example to support 


this theory. 
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3.DEFINITIONS AND NOTATIONS 


A.FREQUENCY PARAMETER 

The frequency parameter is a quantity which indicates 
the position of the driving frequency relative to the 
resonance frequencey,f,,of driven mode in terms of the 04 
of the driven mode.The frequency parameter is defined by 


F172Q4 CE-£4)0/f4 Gr) 


B.STRENGTH PARAMETER 
The investigation of the pressure waveform in the cavity 
required the calculation of the strength parameter from the 


observable quantities.The strength parameter is defined as 


STRPM=M BQ, 929) 
where 
= 2 
Mg * P4 / CA3Cg ) Cai.) 


and P, is the peak amplitude of p,,the pressure distos55tron 
of the driven mode. 
Е-е Ср орсегуасрте i eal culable quantities it is 
reformulated as[5] 
STRPM= V7 VBQ,/(S_ PCy) (3.4) 
where V and a are the rms voltage reading and microphone 
sensitivity respectively of the receiver used to sense the 


standing wave. 


D 








C. e, is defined to indicate the position of M relative 


Bomene classical harmonize Frequencies ,nw,. 
e 





е - (3.5) 
n ПМ1 
and one can relate e, with n luam c2.18) Such as 
F_=20_e Cac) 
n mn 


D. A pictorial representation of Er which will be useful 
throughout the development is given in Fig.l.From now on 
three subscripts will be used for convenience. i.e. 


F becomes F ; 
n Inu 


Qami 
= ( Zn mt 2 
nw 
FIGURE 

f nm? E 

am o O INE 9 i 
nml 
"nml nn 7,22 "1/2" °nml 0.8) 
nmel 


18 





КЖФИРОКЕТИФАТ DEVEEOPMENT 


O OVT IY DESCRIPI Ven 

Assume a perfectly rigid-walled rectangular cavity which 
has one wall perturbed such that the cavity dimensions are 
L [1*€ £Cy,2)] ШЫ and L, as shown in Fig.2 below.Also assume 


the perturbation on the boundary 16 very small compared to 


imemeayity Gimensionss |e2t,,z) | 1. 
z 





FIGURE 2 
The cavity 1s to be excited by a source near the origin 
in such a way that the (N,M,L) mode is driven at a frequency 


close to its resonance frequency. 
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Peele PERTURBED BOUNDARY 
For a rigid-walled rectangular cavity with ideal boundaries 
(€=0),the pressure Po obtained from the linear wave equation 


neh losses 
2 


LI. р0:0 (Ha i) 
MS СС to the following conditions, 
Vp, -n=0 ат = le Cu?) 
X 
ус 
zs 1], 
2 


where n is the local normal to the ideal boundary нее а Еол 
кор Po in terms of Mach number is given bylu] 


Bio 





„My cosk x cosk y cosk,z EOS (NETO ) 
De 7 (4.3) 


-Mg ОМ. БӘ) 


кесе k's are given in Eq.(2.2), and 


o “ek “ex 551 
a у zZ 


If the cavity has perturbed walls,the solution will be 


(4.4) 


em ог а Summation of the Classical linear solution 
for ideal boundaries plus perturbation correction terms 
Eur Eco erregular boundary: 


оле se C4. 5) 


p=pot€p'+€ 
Since the magnitude of the boundary perturbation 1s kept 
small,second and higher terms in € can be considered 


insignificant,so that 


p-pg*€p' (to first order) (uc) 
and satisfy the following conditions, 
P= (4.7) 


29 





Yp.A=0 at E | (4.8) 
LIE. 
z=0,L 
2. 


where N,the local normal to the real surface,is obtained by 


taking the gradient of the equation for the boundary,given 


by [2] 


3= yfx-L Liresty,2]) CEW) 
MWe to the first order in Є, 
h=x-€L Ор) ет, f2) 2 Ce) 
X OY x 02 
КО Мое кс. 105 ле used@im (4.8) the result is 
E 28 OS) ӘР E f2) 4 =Ü s 
ox x 799 dy 92 92 xely li +efon2)] 


Pau iss —expanctomis| tor р evaluated at the real 


boundary L [1*6f(y,z)] produces 


+ ۴ E С(У,2> 
= х 
Pie se cise f0y.29) Pl. Lx ui ail 
A (У, 2. + „к еы (4.13) 
Т 2 ae uuu сє: Fo 7" 


Ae ЛЕШЕ Eq. (A. 6) into Ras of (4.13),taking the partial 
derivative with respect to x on both sides and keeping the 


first-order terms in€ ,yieids E م‎ 
9 


IX IX i Сі РСу,2) +: -- 
OX | аа 
x=sLyLUi i efty:2)] ha. aly OX 


х=, (4.14) 


Taking the partial derivatives with respect to y and z and 


using exactly the same procedure gives 
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DY ОМ ire] xily x-Lx Raby 








- 9% 
x= ee 22 
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-" pe 
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Ex “02 Ox 


dP 
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21 f (Y2) T 
хіх 

















(4.16) 
ET EUCLID C.C) CEI O ana (4.16) into (4.12) and keeping 


ШЕШ (5 Order terms їп ЕЁ results in 


Әр! к b em 27%, FOD 2 , Fz Ifo 
2X 'xzL, "me 93 9) Qz 92 -xzly (, 17) 





The RHS of Eq.(4.17) can be represented as a Fourier 
series in cosines,so that p' can be expressed as a summation 


of normal modes. et 


LT 
=>. Am cos Ty cos 77 cos(wt+0) 
Х. = ыт mz 


or (UE 
co 


- Y > am (0 m,l/wı8) 


= Mey mz0 1:0 








where Ec 
Ano dm JJ [6] dydz (4,193) 
yz 
а 
dmo = on )) [G] cos — y ول‎ dz (4, 19D) 
А іт 
dor = = {J Er dy dz (4.190) 
y 2 
| ЕШ. m+0 
3 ЕЕ |) CAET cos T 2 dy dz ›+ {+0 (4.19d) 
yore 2 
and G is defined as 
2 OP 
= 5 b. (4,19e) 
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Tn oro ПЕК Ее вопегіЬритіоп to p 


from each of the 
terms,it is stated[?]that for a cavity forced by a dynamic 


boundary condition at a boundary 
2 











e "un" p'=0 (4.20) 
and the dynamic boundary condition from the (m,l)th term is 
! s 
22, L(wt +8) 
( mm LIT 
= Aces m е Cae ap) 
Ox TA y = 
then 
| ІТ ¿(wt+09+679mu) 
2 c. 2 As Sue X COSTI cos l— mo 
Б de yx (4.22) 
where 
my 1 20 т 
A if n=1,2,3, B 
and S ml emo Ven soy roc. и 


lS UNIS soltero. Eq. (+.18) gives the complete 


solution. for the first-order perturbation 


p= ыч 


oo 


(4,24) 
— > Sim Or) Spmi (ni m.L /w, 9 + Cam) 

077, 

i 

and combining (4.24) with (4.6) yields the total acoustic 
pressure in the Cavity 





? Z(NM,L/w,8)- (x SL 24, „дл (1) Sam (9 mil /w, 0+ Cane) 


er 
Nee 
Шәміс sos mean a resonanee .wayW 


‚then this term will 
nml 
dominate the summation and all other non-degenerate 


Terms 
сап be omitted.Consequently,the Eq. (4.25) becomes 





n 
P- (mmL/wme)- & Ami 4n(-4) S 


me Cm L/w, 9+ Tee) 
eN x 


СОВ) 
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BERSPECHEIC EXAMPLES 
Шы САЛУ ИЕГЕ STEPPED PERTURBATION 
Assume that the rigid-walled rectangular cavity given in Fig.2 
is perturbed as shown in Fig.3 below,and also assume that the cavity 
is driven in the (0,1,0) mode resulting (0,1,0/w,0) standing wave 


and that the (0,2,0) and (1,0,0) modes are (nearly) degenerate. 





FIGURES 
From Fig.3 the equation for the boundary at L, can be found, 
x etie S [16-10 - 10-0]] (sa. 
By means of Eq.(4.8) € and f(y,z) can be written as 
E МЕЕ (SA. 


Bes] at txt (5А. 


Since the (0,2,0) and (1,0.0) modes are degenerate the emphasis 
of this development will be on these particular modes.The pressure 


i T Diurion ОЁ (0.2.0) mode is 


Р, cos 2T y cos(2wt+ 82 ) 
Ly 
or (5A. 


"P520 Р, (9,2,0/2м/,62) 


where P, is the amplitude of (0.2.0) mode. 


Ро = 


u 


20) 
20 


ц) 








Utilizing the theory developed in the preceeding sections and 
using the equations (4.17) through (4.26), the first-order perturbation 


correction can be found 





Op! E = z) O0 P20 
= BR" LLS =] m (5А.5) 
sime 
[pog 2T y cos(2wt 02) ]| 6(y-U') - &(s- لد‎ 
Ly Ly 
ЕО 0 Сап Бе written as à Fourier series | 
9f. _ 207031, à. cos у cos(2wt +07) (54.6) 
Ox xl, ТЕТЕ 5 а ЕЕ; 


E In GE Ee E (G25) MII CE J vields the Fourier coefficients 


OL PAR TS mm (') - sin( 20.) cos( ML) | 
and 
L [sia (2E L') - sin (22 C) ] 
do q sin ( Ly - ) (5A.8) 
Recalling Eq.(4,23) and (4.24) first order perturbation correction ;p. 
is found as | 
1(2м%-0,% 600) 
T E 
y ET, 
and | , 
c( 2Wt * O2 + C100 
EP = 4T Pz га. L —L— , Quoc sin Ono cos TE е SAID 
BE. (са х 
9 
where 
0 2 2 2 
2W M se 
(22) = (2 Rigg ) ч Ey ) 
Hence ,the total pressure,associated with the angular frequency 2w, 
in the cavity is 
7 
F = Fozo + EP 
- B (0, Z 0/2w, 92) то ш Ша € d Qoo Sin Go 100 c i. 0, 0/2м/, Өз бы) 
(САСТЫ 


2 





The total pressure at the microphone position,x=0 апа у= „is 


| у 
«(2wt4 92) | ((2у++Ө 700) 
Ph. ion 2 eje + 2 €doQ100 Sin Go] € 3 | (IO) 


Define s-[ ] and after a little manipulaticn and use of trigonometric 


identities ,(5A.12) becomes 


EL esr ОСЫ Fon ) Cos(2wtt 92) - (B sin Too ) Е) 
Ger) 


s crpintude of the total pressure in the cavity is 


2. 
P mic. position = P (1+Bcos Two) + ( B Sin Ooo) (5А.1%) 


Earl) iseine corrected vadue ot the amplitude of the second 





harmonic of the driving mode ,obtained by Eq.(2.16). because of the 
boundary irregularity given in Fig.3. 

Now,it is desired to express sin Gr, in terms of the frequency 
parameter of the driving mode,(0,1,0).With the help of Fig.1l.sin@5, 
can be written as 


À 1 
Sin c. Gael oy) 


= f 
% [pSEI- SR) T] 


If f approaches to zero then QG.. approaches toT ,and if f 


100 


approaches to infinity then Q4 1s close to zero. In these same 


00 


limits cos Gio ë 


T 
BOS Ga. = Ji- и, (5А.16) 


goes to -1 and +1 respectively. Hence 
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For 2562 120 ,(5A.15) becomes 
Í. 


sno. = — S 
100 DES = E БАТ” 


pone Eg.(3.l1) and (3; к Fun and e4gg can be written in the form of 


eg el) (5А.18) 
E Oto 
Cinco = fwo - 2Гою CSA IS) 
24010 
Eq. (54.19) can be solved for 1۵0 ЕНЕР іле substituted INO (SA. 17) 


26. 2 £ - fioc s 208 Qioo £-fo10 (1+ 8100) 





AZ 0) 
100 Чоо Foro (1+ C100) 
Use of 1/(1*e,gg)7-1-e49g and little manipulation reveals 
| | 1 
абы Z ——— a nn (5A.21) 


y! + | азе Ba 2 Qoio Coo] (1- tw). 


As a result,equations (5A.15).(5A.16).(5A.21) are the final amplitude 
correction of the second harmonic of the driving mode obtained by Eq.(2.16) . 
A computer program for this was developed by author and is given in 


appendix A. 


x 





B. LINEARLY PERTURBED CAVITY 
Using the same assumptions in section A,assume that the rigid-walled 


rectangular cavity is perturbed linearly as shown in Fig.4 below. 





J 
Ly tomate 
+Ly | 
m 3 
FIGURE 4 
The equation for this perturbation is 
E eres  ( E 
X n = ( EJ (58.1) 
Hence, 
с --2. (5B.2) 
Lx 
апа 
SED 


2 (1-<£ y 
fy2)= (1-52) 
Applying the seme procedure as in section A,the first-order perturbation 


correction and the total acoustic pressure associated with angular frequencv 


2w in that particular cavity can be found 

L J[- 27& sis 2T, cos (2wt «6, ) ] (- | 
«it Ly Ly га 
Clee) 


2T 
OX Lay 


T : 
EY зіп 27 у со5 (204+ Ө:) 
E Ly 
I 
Fq.(5B.4) can be written as a Fourier series and the Fourier coefficient, 


a .is found bv an integration procedure evaluated in the interval 





Ü ta > .Ihe result is 
дш “-+] eos (2-M)T ‚ соз(2+ту — |  — | md? 
Tl (2-m) (2+m) (2-т) (2-т) 
do = 0.0 
SEE) 


а: - 0.0 
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Recalling the Eq,.(4#.24),the first-order perturbation 


correction for (0,2,0) mode is 


mei d і 
f = "^ x Zen, ССОО (1.0,0/2W, 03400 ) (5B.6) 
Nc x 


and the total acoustic pressure associated with angular 


frequency 2w in the cavity becomes 


- (0,2,0/2ш,6,)- “TÊ 69 (2)(4) Sioo (1-0,0/20,02+*000) (58.7) 


12 ZEN! 
I ic 


since agz0.0 

2- (0.2,0/2м,9:) (58.8) 
According to the calculation developed above there is 

no need to make a first-order perturbation correction to the 

E mode іп the cavity Shorn in Fig.!t. As a result,the 

pressure distribution is equal to the second harmonic of the 


driven mode,since a¿=0.0 and this yields p'=0.0 
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БТ ПЕПСЕІ РЕСТОВВБВАТТОМ 
Under the same assumptions made in sections A and B, 


assume that the cavity is perturbed as shown in Fig.5 below. 


y 
Ly Aa 


«ме un өше g 


Ib 


Lx 
FIGURES 


The equation for this perturbation 15 
ші -4/23 
X a = ( Ly 1) | Сею, 


ЕЕЕ EC ЕТЕ 272 20022) сап be written as 


EA 
pe (56-2) 
апа 
Го,2) = 0.0 , Y Ly/2 
2 ac 
p" , Y 3L3/2 


Applying exactly M same procedure followed in section A, 
the total pressure amplitude in the cavity(in first-order 


Peteurbation) is 


P| = P, V(I+ Beosóo ) + (B sin Go). BE 


mic. робат 


where 


B = —— 99 € Quo Sin Groo 


sin Gio 


The theoretical predictions of these specific examples 


апа совСс are given by Eq.(54.21) and (54.16) respectively. 


were examimed with series of experiments developed by (8 | 


The further discussions about these will be given in the next 


soe ron. 
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OTE ULIS 


In this section the theoretical predictions performed 
in sections 5A, 5B and 5C will be compared to the experimental 
results obtained by [8]. 

The information on the empirical losses and resonance 
Frequencies is contained in the Q's and e's.These are the 
values used in the computer program to predict the harmonic 
КИСЕ r TOO the basis Of Eq.(2.16) and is plotted as thin 
ss сируев.ІПһе геспійіс of including the first-order 
perturbation correction are plotted as thick solid curves 
for each specific example.The theoretical curves in figures 
6 through 16 were plotted along with the experimentally- 
measured values for the cavity configurations shown on top 
of each figure.The theoretically-predicted values were 
generated for frequency-parameter intervals of 0.2, and 
the experimentally-measured values were plotted as square- 
blocks.Data were taken, and theoretical predictions made, 
for different strength parameters for the (0,n,0) mode 
associated with different cavity configurations.It is 
important to note that the n=2 distortion peaks when the 
system is driven at this frequency.That is, when the driving 
frequency w is equal to (1/2) w,,there is maximum content 
Sie eines point: where this occurs for each POP curve is 


2 j. 


gn : 
indicated by¿arrow with the label of Fooo At that point 


the value of frequency parameter is 
"02074 %010 Pozo 


(The same thing could be done, of course, for any member 


ou 








of the (0,n,0) family). The arrow labelled as Pr00 Madicates 
the position where the nearly degenarate (1,0,0) mode is 


resonant, and the value of P199 15 


ee eo “100 

The theoretical predictions made in section 5B were 
compared to the experimental results as seen in Fig.7. 
ИЮС Бісв.) 1s compared to Fig.6, the unperturbed cavity, 
it is clearly seen that the theory and experiment are 
excellently in agreement. 

Боп а wedged perturbatien, the theory predicts the 
Beeaueney Gf the second harmonic at which the effect of the 
perturbation occurs as seen in Fig.8.The predicted magnitude 
SE perturbation effect fer this configuration is in good 
agreement with the experiment.The anomolous behavior of the 
posusb»spuonsc 1n Fig.8 is unexplained. 

For stepped perturbation, when the cavity is perturbed, 
but the geometry leads to no predicted correction as seen in 
Ee Eesds EO predicted Correction less than about 0.02 
КЕРЕ СГ less than арвопе 8.05 ас іп Ғіс.1Ш%, then it 
was observed that there was very little or no effect from 
the (1,0,0) mode.Agreement for these cases is good except 
ior the region lying between frequency parameter 4 and 9 
in Fig.9.What happened in that region is also unexplained, 
but it was observed one time only.When the amount of 
perturbation correction is increased the theory predicts 
effects larger than experimentally observed.However, the 


errect of the perturbation appears at the right frequency 


5/2 





parameter as iS seen in FPrg.'s 10,11,13,15 and 16. 

Choosing the shim position, length and magnitude is very 
important as well as is choice of the strength parameter. 
For the shims, M=0.04 and 0.25 inches for stepped and wedged 
perturbations respectively.The effect of strength parameter 
can be seen in Fig.'s 15 and 16.The experimental data associated 
with the third harmonic in Fig.15 were believed to come from 


n on l dıstortıion in the piston motion. 
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FIGURE 16 








Ue CONCLUSION 


Non-linear theory has been applied to standing waves ina 
rigid-walled rectangular cavity with a perturbed boundary 
in order to find one possible mechanism for the excitation 
of a standing wave other than those belonging the family of the 
driven mode.It was observed that such an excitation exists 
if the boundary perturbation and the dimensions of the cavity 
are favorably chosen. 

It appears that the present theoretical model succesfully 
predicts the major features of harmonic content for finite- 
amplitude standing waves in the cavity when the geometry leads 
Benno perturbation correction tEig.? and 9).When the magnitude 
of the perturbation is increased the predicted features were 
larger than experimentally observed.Second-order perturbation 


corrections may be needed to account for these discrepancies. 


5 








The original computer program for Eq.(2.16) was prepared 
by Coppens in 1973, and author made some extentions to that 
eer mM SO that 1t would (1) calculate the perturbation cor- 
Beetion and (2) present the results graphically. The program 
calculates the relative amplitudes and phase angles of sten- 
ding waves in cavities keeping the strength parameter constant 
and changing the frequency parameter to generate response 
curves showing the amplitudes of the nonlinearly excited stan- 
ding waves as function of the frequency parameter ПЕ туссе 
el eS the perturbation correction according to Eq.(4.24) 
and then finds the total relative pressure amplitude using 
EL О ite program alse draws the graph of the relative 
pressure amplitudes of the ideal cavity,total relative pressure 


- 


amplitude of the perturbed cavity and the experimental data 


P 
4 


оп а 3 cycle semilog paper.The Versatec Graphics Plotting 


i-t) 


Manuai f10| was used for the graphical processes on the ІВМ 


да 


360 ог the Randolph Church Computer Center,Naval Postgraduate 





SOME USEFUL INFORMATION ABOUT COMPUTER PROGRAM 


Quantities marked with (*) must be controlled or changed 


hor each run. 


*KON 


*NCUR 


=NDAT 


BUR(I,J) 
*XL 
AY], 
*DELTA 
*STRPM 


*F REQ 


XDAT (JET) 


YDAT( JET) 


The number of iterations throughout the region 
of interest.For this program the iterations are 


performed with 0.2 intervals. 


The number of the experimental curves to be 


drawn + 1 


The number of experimental data in the region 


of interest 

The array that stores the experimental data 
Ilo nm сЕ Л Улгу in the x-direction 
The length of the cavity in the y-direction 
The magnitude of the perturbation 

Strength parameter 


Frequency parameter stored in ATA(I,1) and 
ER TT DUC E TNE maximum value of FREQ 


in the region of interest 
The x coordinate 


Оле ot the eun emos) for ХВАТ 


ATACI,N) ,N#1The array that stores the logarithmic value 


of the pressure amplitudes of the harmonics 


of the driving mode 


ZER(I,N),NZ1The array that stores the linear value of the 


OUI 


SEI 


EG е Оа s Kerr rh ohasmonies of the 


driving mode 


Quality factors of driving mode and harmonics 
orat 


Сыз ЖОО е сыл Ге апсі пагпопісв of it 
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“0100 
*E100 


Е CJIETt1) 


D T O JET+1) 


XDAT(JET+2) 


IDAT JET???) 


HUM 


elo 
*B 


€ 
= 100 
= 0100 singjigg(ag/2) for wedged perturbation 


MINIT сасзеш ев CT 70m0) mode 
e-value of (1,0,0) mode 


Integer value of left-hand corner on the x-axis. 
It must have the same value as the 7th argument 


of subroutine CALL AXIS for x-axis. 


Integer value of left-hand corner on the y-axis. 
It must have the same value as the 7th argument 


of subroutine CALL AXIS for y-axis. 


Increment value of x and y for scalling purposes 


The linear value of the total acoustic pressure 


amplitude associated with angular frequency 2w. 


SinGy),61/4) ur, a, for stepped perturbation 


and ag 15 thewD dE curyer coefficient. 
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